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We report on optical imaging studies of self healing after laser-induced photodamage in disperse 
orange 11 dye doped into PMMA polymer. In particular, the high spacial image of the damage track 
made by a line focus pump laser allows the recovery rates to be measured as a function of burn dose 
using the relationship between transverse distance and pump intensity profile. The time evolution of 
the damaged population results in an intensity-independent time constant of r = 490.2 ±23 mins, in 
agreement with independent measurements of the time evolution of amplified spontaneous emission. 
Also observed is a damage threshold above which the material does not fully recover. 
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INTRODUCTION 

Photodegradation is a processes by which light dam- 
ages a material, from sunlight fading colors on upholstery 
in a car, to catastrophic damage to a material in an opti- 
cal device caused by a laser. Over the past two decades, 
research has focused on understanding laser damage in 
order to find ways to make materials more resistant to 
such damage, [1, 2] and to extend the life of laser compo- 
nents or waveguide polymeric devices. [3-6] 

Recently, using different methods, a variety of dye- 
doped polymers were observed to decay and then re- 
cover from lasing damage when left to rest. Using 
fluorescence measurements, Peng et al reported that 
rhodaminc-doped and pyrromethenc-doped polymer op- 
tical fibers recovered after photodegradation. [7] How- 
ell and Kuzyk [8, 9] used amplified spontaneous emis- 
sion(ASE) to observe decay and recovery of l-amino-2- 
methylanthraquinonc (disperse orange 11, or DOll) dye 
doped poly (methyl methacrylate) (PMMA). ASE was 
found to recover fully even when the damage of the dye- 
doped polymer was nearly 100%. However, no recovery 
was observed in liquid solution of DOll in MMA - the 
liquid monomer used to make PMMA polymer. 

These results suggest that the host polymer plays an 
important role in self healing. While orientational hole 
burning due to molecular hole burning is a common phe- 
nomenon observed in such dye doped polymer systems 
[10], no dichroism is observed at the higher pump in- 
tensities used by Embaye and coworkers, which rules out 
molecular reorientation as the underlying mechanism. [11] 
We also find in our present imaging studies that there is 
no dichroism associated with the burn track. Embaye 
also complemented ASE measurements of DOll-doped 
PMMA thin films with linear absorption spectroscopy, 
which showed that during damage, the DOll molecule 
is transformed into a quasi stable damaged state, as sug- 
gested by an isosbestic point. Zhu et al. [12, 13] used 
two-photon fluorescence of AF455-dopcd PMMA poly- 
mer to measure decay and recovery, and confirmed that 



this system also recovers fully after photodegradation. 

The present work seeks to study the evolution of self- 
healing of DOll dye-doped PMMA thin films using a 
change in the absorbance on the damage track as a proxy 
for the damaged population. Since the excitation line has 
a gaussian profile, such a beam will create a burn mark 
with a gaussian damage profile. This allows the study 
of recovery rates as a function of the degree of initial 
damage by imaging the back-illuminated burn mark of 
the sample using blue light. From the time dependance 
of the intensity at each pixel, the recovery rate at the 
corresponding position of the sample is determined. By 
exploring recovery as a function of dose, we can investi- 
gate whether or not the recovery process is independent 
of the degree of damage; or, if there is a correlation be- 
tween the amount of damage and the rate at which the 
system recovers. Such information is an important input 
into the determination of the underlying mechanisms. 



EXPERIMENTAL SETUP 

Photodegradation and recovery is studied in 
DOll/PMMA using optical imaging. The appara- 
tus consists of a computer controlled Nikon Coolpix 
990 digital camera with a custom microscope objective 
mount, a custom sample holder mounted on a three 
axis translation stage, and a blue LED (Radio Shack 
276-0316) peaked at 463nm for controlled sample 
illumination. The spectrum is shown in Figure 1. An in- 
dependent UV-Vis absorbance measurement shows that 
the largest absorption change due to photodegradation is 
observed in blue part of the spectrum. So, to maximize 
the observed change in transmittance, we probed the 
sample with a blue illuminating LED. The LED is con- 
nected to a DC power supply via a 97.7 ± 4.3fc57 resistor 
and the applied voltage is maintained at 4.49 ± 0.01U. 

The camera was set to use the sRGB IEC61966-2.1 
color scheme, with an FN = 4, and an exposure time of 
1/60 s. A 20x objective was attached to the camera and 
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FIG. 1. The LED emission spectrum. It is sharply 
peaked(cr = 14.37 ± 0.09 nm) at A = 463 nm. 



the camera's maximum optical zoom was set to 3x. By 
imaging a fixed-width line, it was determined that 1 pixel 
= 1.713 ± 0.062/im. For simplicity pixels are used as the 
unit of length in this paper. 

The sample used in this study is a 9g/l thin film of 
Disperse Orange 11 dye doped in PMMA polymer that 
is sandwiched between two microscope slides. First, dye 
is dissolved in MMA and placed in an oven at 95°C for 
2 hours to form a highly viscous, semi-polymerized so- 
lution. After cooling to room temperature, a few drops 
of solution are placed on a microscope slide and spread 
across the slide using another glass slide to form a thin 
film - a process known as doctor blading. Two such slides 
were prepared and placed back into the oven for a few 
more hours. After being fully polymerized, the slides are 
pressed together - polymer to polymer - and placed in a 
custom thermal press to form a uniform film. The press 
temperature is maintained at 150°C and constant uniax- 
ial stress for two hours resulting in uniform thin films. 
The sample used in the present study is 150/im thick. 

The thin film is photodegraded with a frequency dou- 
bled Nd-YAG laser with the following specifications: 
532nm wavelength, 10ns pulse width and operating at 
10Hz frequency. The beam is focused to a line using a 
cylindrical lens so that it has a Gaussian profile across 
the line of Gaussian width a = 51.67±0.03/zm. The sam- 
ple was photodegraded by exposing it to the pump beam 
for 31.4 mins with a per pulse energy of 90//J. Details 
of the experimental setup for burning the sample can be 
found in the literature. [11] 



METHOD 

After the burn was complete the sample was immedi- 
ately mounted in the imaging apparatus where it was 
aligned and brought into focus for the camera. Im- 
ages were then taken using exponential time intervals 
beginning with two-minute intervals, followed by fifteen- 
, thirty-, and sixty-minute intervals. The images were 
then imported into Igor Pro where the RGB line profiles 
were analyzed at several cross sections with an average 
burn width of 20 pixels. Figure 2 shows an image of a 
burn spot. The burn was found to have a Gaussian pro- 
file and was most pronounced in the blue channel, which 
agrees with the spectral region to be most affected in in- 
dependent linear absorption measurements. Using just 
the blue channel, we fit the Gaussian profile for each im- 
age as shown in Figure 3 for a few representative times. 
Using the Gaussian fits we determined the value of the 
new parameter j3(x), which represents the normalized in- 
crease in transmission at position x along the burn, and 
defined by, 

m - (i) 

where Bq is the baseline transmitted intensity, and B(x) 
is the blue channel intensity at position x along the burn. 
We define the normalized increase in transmission at the 
center of the burn as 



where A is the amplitude of the Gaussian burn profile. 
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FIG. 2. An example of a burn line used in this study(edited 
for clarity). Cross sections were taken perpendicular to the 
burn line with a width of 20 pixels. The dark spots are the 
scattering patterns from micro bubbles in the polymer. 



THEORY 

In this section, we relate the imaged intensity to the 
population of damaged molecules. To determine /3 from 
the image, it must be recognized that the blue channel 
intensity, B, is a convolution processed in the camera, 
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FIG. 3. Gaussian burn profile of the transmittance as a func- 
tion of position at several times. The data has been statisti- 
cally smoothed for ease of viewing. 
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FIG. 4. Four images of burn lines over time(contrast and 
color edited for clarity). 1: t = mins, 2: t = 600 mins, 
3:t — 2739 mins, 4: t = 5439 mins. Note the visible dimming 
of the transmittance profile as a function of time 



given by, 

B{t,x) = J S{X)I (X)T(X;t 7 x)dX, (3) 



where S(X) is the sensitivity of the camera's CCD in the 
range X m in < A < X max . Io(X) is the LED intensity 
spectrum and T(A; t, x) is the transmission coefficient of 
the thin film at time t, and position x. 
From the Beer-Lambert Law we have, 

T(X;t,x) = e-^M*-*)' (4) 

where er(A) is the absorption cross section, n(t, x) is the 
number density of the absorbers, and / is the thickness 
of the thin film. 

We base our models on observations from photodegra- 
dation and recovery studies of ASE in DOll/PMMA that 
include a long-lived damaged population and an undam- 
aged species that can become damaged through photo- 
excitation. [11] The fractional undamaged population, n u , 
is given by 

n u (t;x) = 1 - n (x)e~ t/T (5) 
and the fractional damaged population, rid, is 

n d (t;x) =n (.T)e- t/T (6) 

where no(x) is the fractional initial damaged population 
at position x along the burn, and 1 /r is the recovery rate. 

Recognizing that the damaged and undamaged popu- 
lations have different absorption cross sections (cr'(A) and 
cr(A) respectively) we can substitute Equations 6 and 5 
into Equation 4, 

T(X;t,x) = e - (T ( x ) l ( 1 - n o{ x ) e ~ t/T )~<?'W ln a(x)e~ t/T ^ ^j) 
and then, rearranging slightly, we get 

T(X;t,x) = T {X)er ArT ^ no ^ le ~ t/ \ (8) 

where Act(A) = cr'(A) - a(X) and T (A) = e^^ 1 . Sub- 
stituting Equation 8 into Equation 3 we find 

B{t;x) = / S(A)/ (A)To(A)e- ACT ( A )"°^ e " /T dA. 

(9) 

If Jo (A) is sufficiently narrow and centered at Ao (as is 
the case for this experiment) we can approximate it as a 
delta function and Equation 9 becomes 

B(t; x) = S(Xo)Io(Xo)To(Xo)e- A ^ n °^ le ~ t/T . (10) 

Substituting Equation 10 into Equation 1, we get the 
transmittance as a function of position and time, 

P{t;x) = Ae- a{x)e ~ t/T 7 (11) 

Where A is related to the sensitivity and transmitted 
intensity, a(x) = Aa(Xo)no{x)l, which is related to the 
change in absorbance, and 1/r is the recovery rate. 
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Assuming that the degree of damage is proportional 
to the beam intensity, which has a Gaussian profile, a(x) 
will be of the form 



a(x) = a(xo)e 



-C- 



•-f 



(12) 



where xq is position of the center of the burn, a(xo) the 
change in transmittance at the burn center, and a is the 
beam width. 

Because the change in absorption is small relative to 
the background intensity, the intensity difference deter- 
mined directly from the raw image is noisy. The intensity 
profile for further analysis is determined by the values re- 
turned by the fitting function. Thus, the intensity at each 
pixel position is determined by first fitting the image data 
to a Gaussian, and then using the fit parameters xo, A 
and a to determine j3(t;x). 



DATA AND RESULTS 

We first studied the recovery process at the center 
of the burn, /3(f;xo), which is defined as p in Equa- 
tion 2. Using Equation 11, we performed multiple fits 
with varying initial parameters and found r = 490 ± 23 
mins. These results arc consistent with literature val- 
ues of t — 476 ± 14 mins, which were determined from 
decay and recovery studies of Amplified Spontaneous 
Emission(ASE).[ll] Our only assumption in making this 
comparison is that the ASE signal is a measure of the 
undamaged population, or Iase oc n u (t). 
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FIG. 5. The degree of damage as determined from the peak 
burn dose at the beam center, p, as a function of time (points) . 
The fit (curve) yields r = 490 ± 23 mins. 

We can explain the small discrepancy in time constants 
by noting that the ASE intensity, Iase, is the total inten- 
sity that originates from all parts of the excited sample. 



Thus, the ASE intensity is spatially integrated over the 
full burn profile, or 



(13) 



lASE{t) oc / n u (t,x)dx 



If we consider the area under the Gaussian burn as a 
function of time, as shown in Figure 6, and fit it to an 
exponential (which is the equivalent procedure used for 
the ASE measurement), we find r = 472.3 ±8.8 mins. 
This is in agreement with the ASE measurements. 

The time constant, r, which is determined from the 
time-dependence of the increased transmittance at the 
center of the burn mark profile, p, gives the most precise 
measurement of the recovery time due to the fact that 
the change in transmittance there is the largest. 
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FIG. 6. The transmitted intensity averaged over the trans- 
verse profile of the Gaussian burn as a function of time 
(points). The exponential fit (curve) yields a time constant 
of r = 472.3 ± 8.8 mins. 

Next we determine the recovery rates as a function of 
distance from the peak to determine the recovery rates as 
a function of dose. Since we postulate that the recovery 
rate is independent of dose, the time constants should be 
independent of a(x) since it is proportional to the dose. 

Given that the image is noisy, especially in regions 
where the dose is small, we use fits to the burn profile 
to determine the transmittance intensity at each pixel, 
as follows. At each time, we fit the data to a Gaussian 
to get the amplitude and width. Then using the fitted 
amplitude and width, we plotted the transmittance at a 
given a; as a function of time, see Figure 8. By fitting this 
data to our model for the transmittance, Equation 11, we 
find that a(x) is indeed Gaussian (see Figure 7) with a 
width of a = 28.8 ±1.9 pixels or a = 49.3±3.7/im, which 
is within one standard deviation of a direct measurement 
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tion of time for several distances from the center of the 
burn line as quantified by the pixel number. The inten- 
sity at each pixel is determined using the fit parameters, 
as described just after Equation 12. The solid curves 
are the associated fits to the data of the master function 
/3(t;x). Since that data is noisy, only three sets are dis- 
played. Where data is not shown, the fit functions are 
shown. Fitting f3(t;x) to each of the decay curves deter- 
mines the recovery time constant. The average time con- 
stant from these curves is found to be r = 483 ± 122mins, 
which is within experimental uncertainties of the value 
found for p(t). The large standard deviation is due to 
low signal-to-noise for the lowest intensity curves. 



CONCLUSIONS 



FIG. 7. Using fit parameters we found /3(t; x), which we then 
fit to find ce(x), which is a Gaussian across the burn with a 
width of a = 28.8 ± 1.9 pixels, or a = 49.3 ± 3.7/im. 

of beam at the line focus using a beam profiler. Using 
the fit parameters, we determine the degree of recovery 
as a function of the distance from the center. We can 
thus determine the degree of population recovery at each 
point on the burn profile. 
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FIG. 8. The image intensity as a function of time for several 
pixels (points) and fits to the function /3 (curves). To improve 
plot's clarity, only three raw data sets are shown. Fitted 
curves for the omitted data are plotted. 

Figure 8 shows the change in transmittancc as a func- 



Using optical imaging studies, we characterized the 
recovery of a 9g/l DOll/PMMA as a function of time 
and burn intensity using position from the burn center 
as a proxy for intensity. Across most of the Gaussian 
burn, the damaged populations recovered at the same 
rate, which shows that the rate of recovery is indepen- 
dent of the dose, supporting the argument in our previous 
work that the recovery process is independent of pump 
intensity. 

Further variables to be considered in future research 
are the effects of dye concentration, temperature, and 
polymer composition on the photodegradation and recov- 
ery rates. If the recovery mechanism depends on interac- 
tions between the dye molecules, as suggested in previous 
work, [8, 9, 11, 14] then the recovery times should depend 
strongly on concentration. If the polymer plays a role 
in recovery, changing the polymer or the temperature 
(which softens the polymer) should change the recovery 
rate. 

The observation that ASE appears to fully recover over 
time still needs to be reconciled with the fact that our 
data suggests that the optically visible damage does not. 
Some of the optical imaging data suggests that there is 
a damage threshold above which a material will not fully 
recover. This effect is observed in the tails of the Gaus- 
sian burn mark where the burn fully recovers, but near 
the center, it does not. The level of noise in our current 
experiment makes it impossible to determine whether 
such a full recovery threshold is a property of the mate- 
rial. A more sensitive CCD camera with a higher signal 
to noise ratio should be able to resolve this problem. 
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